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We develop a new variational scheme to approximate the position dependent spatial 

S" 

probability distribution of a zero dimensional manifold in a random medium. This celebrated 

c 

'toy-model' is associated via a mapping with directed polymers in 1+1 dimension, and also 
describes features of the commensurate-incommensurate phase transition. It consists of a 
pointlike 'interface' in one dimension subject to a combination of a harmonic potential plus a 
random potential with long range spatial correlations. The variational approach we develop 
gives far better results for the tail of the spatial distribution than the hamiltonian version, 
developed by Mezard and Parisi, as compared with numerical simulations for a range of 
temperatures. This is because the variational parameters are determined as functions of 
position. The replica method is utilized, and solutions for the variational parameters are 
presented. In this paper we limit ourselves to the replica symmetric solution. 

PACS: 05.40.+j, 05.20.-y, 75.10.Nr, 02.50. +s 
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I. INTRODUCTION 



Recently a lot of attention has been devoted to the behavior of manifolds in random 
media. P|-f2T[| This is partially due to their connection with vortex line pinning in high T c 
superconductors, [fTjffi] but also because of the intrinsic interest in the behavior of interfaces 
between two coexisting phases of a disordered system, like magnets subject to random fields 
or random impurities. In addition mappings are known to exist between one dimensional 
manifolds like directed polymers in disordered media and growth problems in the presence 



of random noise described by the KPZ equation. [22 



The variational method p.[T|.|T6|] appears to be an important tool in approximating many 
properties of the system like the calculation of the roughness exponent of a wandering man- 
ifold in a disordered medium. Recently we have conducted a numerical study of the spatial 
probability distribution of directed polymers in 1+1 dimensions in the presence of quenched 
disorder. [[HJ Directed polymers refer to an interface in two dimensions with no overhangs. 
In addition to the width of the interface, it is important to know what is the probability 
that the interface would wander a certain distance in the transverse direction. We have 
then exploited a mapping of this system into that of a pointlike interface in one dimension 
subjected to a combination of a harmonic restoring force and a random potential with long 
ranged correlations. 

Using that mapping combined with a use of the variational approximation, we were 
able to derive many properties of the spatial probability distribution (PD) of the directed 
polymers. |T?J In particular the function f(a) associated with the decay of the probability 
for different realizations of the disorder has been derived using the replica method. Replica 
symmetry breaking solutions have been used (at low temperatures). The mapping from 
directed polymers to the so called toy-model is discussed in detail in ref. [TO based on 
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previous work by Parisi || and Bouchaud and Orland |J . In Particular the spatial probability 
distribution of directed polymers at large times relates to the spatial probability distribution 
of the toy-model at low temperature. Villain et al. |l[ emphasize directly the relation between 
the toy-model and the commensurate-incommensurate transition in systems with quenched 
impurities. |L| 

Using the replica method the toy-model can be mapped into an n-body problem were the 
interaction between different particles (replicas) is given by the correlation of the original 
random potential. The variational method utilizes an effective quadratic hamiltonian whose 
parameters are chosen by the stationarity method to give a good approximation for the free 
energy of the n-body system in the limit n — > 0. [§j The results it provides for the tail of the 
spatial distribution are not very accurate as will be shown below (Section 4.). In this paper 
we describe a way to directly apply the variational approach to the probability distribution 
rather than to the hamiltonian or the free energy. As we shall see it gives a much better 
approximation for the behavior of the tail of the distribution. For high temperatures the 
results appear almost exact without the need to break replica symmetry. 

II. THE VARIATIONAL METHOD APPLIED TO THE PROBABILITY 

DISTRIBUTION. 

The toy model |l|J^, p^ , p^1 involves a classical particle in a one-dimensional potential 
consisting of a harmonic part and a spatially correlated random part. It is the simple 
prototype of an interface-a point in one dimension- which is localized by a harmonic restoring 
force, and pinned by quenched impurities the distribution of which may be correlated over 
large separations. Its hamiltonian is given by: 

H{u) = ^ 2 + V(cu), (2.1) 



where uj denotes the position of the particle and V(u) is a random potential with a gaussian 
distribution characterized by: 

(V(u>)) = 0, (V(u)V(u')) = o7T~~ \ I" " ^f- 27 + Const. = -f((u - a/) 2 ), (2.2) 

2(1 -7) 

where the brackets indicate averaging over the random potential. In the sequel without loss 
of generality we will take the constant on the r.h.s. of eq. ([2.2|) to be zero. We will make use 
of the replica method to replicate the partition function 

Z = J du exp( -pH ) (2.3) 

and average it over the random potential to express it in terms of an effective n-body hamil- 
tonian: 



(Z n ) = fdojf- du n exp(-[3H) 



1 n 1 

n = + / ( K - " 6 ) 2 ) • (2-4) 

/ 0=1 A ab 

For later use let us also define an n-body quadratic hamiltonian: 

1 n 1 

h = _ ^^l^ob u a uJb (2.5) 

a=l ab 

where the matrix elements of a aa are free parameters at this point. 

We consider cases in which the potential has long-range correlations, (i.e. 7 < 1), and in 
particular we are interested in the case 7 = | because of the above mentioned mapping from 
directed polymers in 1 + 1 dimensions. Let us review briefly the variational scheme used by 
Mezard and Parisi. We start with the well known inequality 

(e A ) > e {A) (2.6) 

The Mezard-Parisi variational approximation || can be obtained from this inequality, by 
choosing 



A{u)x, ■ ■ ■ ,uj n ) = H{ux, ■■■,u n )- h(u!, ■ ■ -,uj n ) (2.7) 
And by denning the average as 

1 r exp ( —pri) 

with 



Tr 



/+oo 
duj\ ■ ■ ■ du n (2.9) 
-oo 



This yields 



nf = -^lnTrexp ( -f3H ) < (H - h ) - - Tr In exp {-(3h ), (2.10) 

where f is the free energy. The variational free energy is given by the right hand side of 
eq. ( |2.1U| ) at the point of stationarity. This procedure provides equations for the variational 



parameters which are the elements of the matrix a^. MP |L3| have solved these equations 
and found a replica symmetric (RS) solution valid for high temperatures and a solution 
which breaks replica symmetry (RSB) at low temperatures. 

Recall that the limit n — > has to be taken when making use of the replica method. An 



important point to bear in mind is that the inequality (2J3) holds only when n is a positive 
integer, provided the parametrization of a is such that Tr exp(—(3h) < oo. This is because 
only in this case the integration measure 

[ Trexp(-/5/i)]- 1 exp (-f3h) du x ---du n , (2.11) 

is a real positive measure on R n , for which the proof of the inequality ( p.6|) holds (see e.g. 
P3[ p. 61). In the limit n — > the inequality can change sign and does not hold in general, 
but one still expects the the stationary point of the r.h.s. of equation ( |2. 1Q ) to give a good 
approximation to the l.h.s. i.e. to the exact free energy. Support for this contention comes 



also from the fact that the variational method becomes exact for a manifold embedded in d 
spatial dimensions when d becomes infinite, and one can systematically improve it by a 1/d 
expansion. P, |II||i7l , |T9 



Let us define 



Tr = J duj\ ■ ■ ■ duj n 8{uo\ — uo) 



(2.12) 



The function 



P h (to) = Tr exp {-j3h ) 



(2.13) 



(with the parameters a a p determined by the stationarity conditions which were discussed 
above) constitutes an approximation to the exact spatial probability distribution, averaged 
over the random realizations of the potential, and given by the formula 

/ exp(-/3#H) \ 



Ph oo) 



Tr exp (—f3H. 



(2.14) 



\Jdaexp(-(3H(cx)) / 

with the limit n — > to be understood. Eq. (|2.13|) has been evaluated by us in a previous 
publication |Tj| and compared with numerical simulations both using the RS and RSB 



solutions for the variational parameters. The result for Ph is [[L6 

'-(3lu 2 \ 



Ph 00 



g 

2nG 



1/2 



exp 



11 



2Gn 



(2.15) 



with 



G 



ab 



(2.16) 



Substituting the RS and RSB solutions for G\\ obtained in ref. ( |R|), we find |T6 



Ph ( oo 



Phi oo 



1/2 



i 27r(l + 7 - 1 T-( 1 +^)), 

1/2 . 

! exp 



27r(l + 7 ), 



exp 



2(1 + 7) 



pfico 2 



2(1 + 7- 1 r-M) 



LO 



T > 1 



T < 1 



(2.17) 



with the 'reduced' temperature T defined as 

f = r ^^im^ 9 y*. (2 , 8) 

In the sequel we refer to the expressions given in eq. fl2.17| ) as "the hamiltonian variational 
approximation" . 

I proceed to derive a new variational scheme which is more appropriate for approximating 
the position-dependent spatial probability distribution. The idea is to use the inequality 
Q2.6|) , but substitute Tr from eq. ( |2.12| ) for the trace in eq. ( |2.8| ). This means, defining 

«A»(u,) = Tr[A{u^.^ n ) exp(-/^H (jug) 

Trexp ( —f3h ) 

The notation (()) is used to distinguish the ui- dependent average from the averaged defined 
in eq.(O). Using A = H - h I find 



P H (uo) = Trexp ( -(3H) > exp (-/3((H - h))(w) ) x Trexp (-(3h) = P v (u) (2.20) 
Again, for positive integers n, and Trexp ( —f3h) < oo, the measure 

[ Tr exp(— f3h)]~ l exp (— f3h) 8{u\ — u>) dui ■ • • du n , (2-21) 



is a positive measure on R n for any value of to and the inequality (|2.6Q holds. For n — > 
this is no longer true, but the stationary point of the r.h.s. of eq. ( |2.2U| ), denoted by 
P v {uj)i with respect to the parameters a a b is expected to yield a good approximation to 
the spatial probability distribution by virtue of analytic continuation. This will be checked 
in Sec. 4 in comparison with numerical simulations. Such a variational determination of 
P v {uj) is obtained for each value of u separately and thus it is expected to provide a better 
approximation than by just using the o a b obtained globally from extremizing the free energy 
(see eq. ( p,10| )), and using eq. (|2.13| ). One should note that it does not follow from our 
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discussion that P v (u>) is properly normalized. Provided it is a good approximation to the 
exact probability distribution, its normalization will be close to one, and correcting for the 
exact normalization will have practically no effect on the behavior in the tail region were the 
probability is very small and where the variational calculation is most important. 

Our next task is to evaluate the quantity ((TC—h ))(uj) needed in order to calculate P v {u), 
see eq. ( 2.20|) . If we use the integral representation for the Dirac 5-function in eq. (|2.12| ): 

dk 



8{u\ — id) 



iku! — ikui 



2tt 



(2.22) 



we can write 



Tr[{H - h) e 



-ph 



dk „■ 



iktu 



2vr 



d&i ■ ■ ■ duo n 



(^J2f(( UJ a-^b) 2 ) + ^J2 a ^a^bj X exp (~^2(G ^abUaUb-ikUij . (2.23) 
V ab ab / \ ab / 

In the Appendix we show how the integrals can be evaluated for general n. The end result 
is: 



{(H-h))(u> 

P (G 
2 



7T 



1/2 f°° 

[Z ab )~ 1 I dp exp(-p 2 

J — oo 



^ + pV2Y ab )\ 



2Z, 



ab 



x/ f -X ah tf 



2j3 + 2(3 V 2/3G n 



G 



ii 



(2.24) 



with 



X, 



ab 



G aa + G 



bb 



2G 



Y a b — (Gi a — Gib) I {X a b 



ab 



,1/2 



J ab 



Y 2 

1 ab 



(2.25) 



For the special case of interest 7 = 1/2 (see eq. (|2~2j ) for /) we find: 
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W -h)){u) = 
~ (3Y ab xH 2 cu 



9^2 

ab 



2G 11 



erf 



1/2 



2/3 + 2/3 ^ ° 



2 Z a },G u 
1 



(3X ab Z, 



ab 



1/2 



2Z ab Gn 



2/3G 



ii 



V 27rG n J 



exp - 



(2.26) 



where erf is the usual error function. Eqs. ( 2.24 ) and Q2.26 ) together with the expression 
( p.2C| ) for P v (uo) constitute the main results of this section. In the rest of the paper we 
consider only the important case of 7 = 1/2. 



III. THE REPLICA SYMMETRIC SOLUTION 

In this section we consider the replica symmetric solution to the variational stationarity 
equations. In this case we search for a solution for which all the off diagonal elements of the 
matrix a ab are equal as well as all the diagonal elements: 

a ab = a a^b 

<?aa = o~d a = 1 • • -n (3.1) 

we also define 

lii = fi + a d + (n - 1)0- (3.2) 

In terms of these parameters one finds 

G aa = — ( 1 + — ] ; G ab = ^ , a ^ b (3.3) 
ft \ H) K 

Thus we have two variational parameters a and fix as contrasted with the hamiltonian varia- 
tional approach were there is only one variable because in that case translational invariance 
dictates fii = fi. In addition, the variational parameters are of course dependent on uj in the 
present case. We define for convenience: 

9 



T = G n + G 12 . (3.4) 

Using the result of the previous section, eq. (|2.26|) , we can express the probability distribution 
P v (uj) (see also eqs.( [2.20|) and (|2.13| ) ) as 



PM = (^) exp {-^ exp ( g(u) ) 

( \ ^ A( g V /2 V J 2/?r \ 1/2 ( ^ 

g{uj) = _ erf - ^ w + exp 



fiiGu J yn^GnJ \ 2fi 1 TG n J 



where the limit n — > has been taken. The stationary points of eq. ( fj.5|) in the a — fi\ plane 
as functions of to, have been obtained numerically for various values of the parameters of 
the model (P,/j,,g), and the results have been used back in eq. (|3.5|) to evaluate P v (u>). The 
results will be summarized below. But first let us examine the simpler behavior in the tail 
of the probability distribution, that can be investigated analytically. 

For large values of u, the error- function in eq.( |3.5| ) can be approximated by sign(tu) and 
the exponential term in the expression for g(tv) is negligible. If we further define the tail 
region as the region for which 

\u\»/3g/fji, (3.6) 
we find that the extremum of P v (u) is obtained for 

Hi n\u\-P9 

^ « [fin + 2p 3 g 2 - 2gi3(i3 4 g 2 + P^) l/2 }/n . (3.7) 

Plugging these expressions back into the expression for P v (u) one finds that the behavior of 
the spatial probability distribution in the tail region defined by eq. ( |3.6j ) is 
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P v (u) « exp(-///3o;72 + #/3 2 |u;| + C) 

2/i Vtt^J 2/ii 2 2 I 2vr i ' v 



where /ii is given by eq. (|3.71 ) . 

The extremum in the tail turns out to be a minimum in the a — [i\ plane, as opposed 
to the situation for n > 1 when it is a maximum as can be seen from eq.( 2.20| ). (On the 



contrary, for very small values of uj we have found numerically, that the stationarity point is 
actually a saddle point.) The fact that the extremum in the tail region is a minimum seems 
to suggest that the approximate expression eq. (|3.8|) constitutes an upper bound to the exact 
behavior. The behavior in the tail given by eq. (|3.8| ) should be compared with the behavior 
obtained from the conventional variational approximation given by substituting 7 = 1/2 in 
eq. ( |2.17| ). Those formulas give a higher value for the tail than that predicted by eq.(|3.8|). 



(This is true even when T < 1 with the RSB solution). More on this in the discussion below. 

IV. COMPARISON WITH NUMERICAL SIMULATIONS 

We proceed to a numerical comparison between simulations and the results of the various 
variational approximations. Again, we limit the discussion to the case 7 = 1/2. We have 
studied numerically a lattice version of the toy model. A suitable interval of the particle's 
position uj ( — 12.5/a//3 < uj < 12.5/yT? ) is divided into 5,000 lattice sites. For a given 
realization, the algorithm generates an independently distributed gaussian random number 
for each site rf, it then generates Vj, the correlated random potential at site j, by summing 
the random numbers in the following way: 

Vj oc sign(«- j)n. (4.1) 

i 

The quadratic term of the hamiltonian is added to this random potential, and then the 
partition function and probability distribution are calculated. 
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We consider three sets of values for the parameters (/3, g, fi) : (0.2, 2 v /7r, 1), (1.0, 2y/7r, 1), 
(10.0, 2.2, 4.6). The corresponding values of the reduced temperature T for these three cases 
are 5, 1 and 0.23 respectively. The data for T = 5 is plotted in Fig.l (solid curve). It was 
obtained from numerical simulations with 50000 realizations of the disorder. The dashed line 



represents the probability P^ given in eq. ( |2.17| ) and derived from the hamiltonian variational 



method. The diamonds represent the approximation developed in the last section. We see 
that it gives a perfect fit to the data for this value of T. We have used Mathematica to find 
the stationary point of eq. (|3.5| ) in the — o plane and then substituted these values back 
into eq.( |3.5|) . The asymptotic formula eq. ( |3.8|) predicts 

P v (u) « exp(-0.1c; 2 + 0.142cj - 2.1) (4.2) 

where we used \i\ = 0.7 from eq. Q3.7|). This behavior is also indistinguishable from the data 
in the range 10 < uo < 28. 

Let us proceed to the case T — 1. The data is plotted in Fig. 2. In this case the noise 
associated with the random potential is apparent. In this case we found it necessary to 
average over 500000 realizations. The data for 50000 realizations is also shown in lighter 
dots. In the tail region it is apparent that the curve corresponding to the lower number of 
realizations has a lower value, since the average is increased by relatively rare events. The 
dashed curve is the single parameter hamiltonian variational fit. The new variational method 
results are represented by the diamonds. The asymptotic formulas predicts \i\ = 0.06 and 

P v (uj) « exp(-0.5cu 2 + 3.545cj - 17). (4.3) 

It is represented by the solid line in the figure. The diamonds do not lie exactly on this 
line, because uj is still not large enough in this range. Again we see that the variational 
approximation gives an excellent fit to the data. 
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The final example is for T = 0.23 using 10 6 realizations of the disorder. The data is 
depicted in Fig. 3. The dashed line is the result of the hamiltonian RS variation. The dot- 
dashed line is the result of the RSB hamiltonian variation which is the appropriate solution 
for T < 1. The diamonds again represent our new variational method results. The range of 
uj simulated is lower than the onset of the tail region given by eq. (|3.6| ) which in the present 



case starts above uj ~ 5, so a comparison with the asymptotic formula eq.( |3.8| ) is not shown. 

We see that the data in this case falls below the result of our variational method, which 
is nonetheless better than the hamiltonian RSB variational approximation. Two possibilities 
come to mind to explain this discrepancy. 

1. It seems quite possible that because the relative strength of the random part of the 
potential is much larger when T is small, as compared to the harmonic part, one needs more 
realizations to accumulate enough statistics for the averaged probability distribution. For 
T = 5 even 5000 realizations already gave us good results. In Fig.l we show the results 
for 50000 realizations, but these are practically indistinguishable from the average of 5000 
realizations. For T = 1, 50000 realizations are not sufficient and we had to average at 
least 500000 realizations to accumulate enough statistics. This is because we have to collect 
enough rare events which contribute to the distinction between the average and typical values 



of the spatial probability distribution. |T5|JT6|| It is quite plausible that to get better results for 
T = 0.23 one has to go to a higher number of realizations. Averaging over more realizations 
(see e.g Fig. 2 for the distinction between 50000 and 500000 realizations) may narrow the 
gap between the data and the variational approximation. 

2. Another strong possibility is that like the hamiltonian variational case it is necessary 
for small T to look for a solution with replica symmetry breaking when using the current 
variational scheme. Such a solution may display a different asymptotic behavior for small T 



than the one given by eq. fl3.8| ). It may also yield a lower value in the pretail region (which is 
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the range of values depicted in Fig. 3). In order to find such a solution one has to go back to 
the full expressions derived in Sec. 2 without making the simplifying assumptions of replica 
symmetry made in Sec. 3, and one has to extremize the probability distributions under the 
more general conditions. The task of finding of a RSB solution is left for future research. 

V. DISCUSSION 

In this paper we have developed a new variational scheme for approximating the spatial 
probability distribution. The replica symmetric solution gives an excellent approximation to 
the numerical data for high temperatures T > 1. For the tail of the spatial probability it 
predicts a gaussian decay with a linear exponential correction. At lower temperature, the fit 
in tail region is not that good, either because the numerical data is insufficient and one needs 
to average over a higher number of realizations, or a solution with replica symmetry breaking 
is needed (or both). Of course since we are dealing after all with an approximation, we are 
never guaranteed a perfect fit. For all temperatures tested the present method gives a much 
better fit to the probability distribution than using the hamiltonian variational method, both 
with or without RSB. 

When the data for T = 0.23 is displayed in a log-log plot, [IIJ a crossover is observed 
from a behavior P ~ exp(— const u 2 ) to a behavior like P ~ exp(— const u 3 ). An empirical 
fit to the data of the form 

P(u) ~ exp(-3.5c; 3 + 20) (5.1) 

is depicted as open triangles in Fig. 3 for the range 6 < uj 2 < 16. Since in this case we 
are not really observing the tail region (according to the definition given in eq. ( |3.6| ) ), this 
behavior may be a transient or intermediate behavior. We should also be cautious about 
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the data in this region, because as mentioned in the last section it is possible that more 
realizations are needed. 

On the other hand, if we do take the cubic u dependence of the log of the probability 
distribution seriously, this immediately rings a bell because of some work done by Villain 
et al. [|l| on the toy model. What they actually showed was that the probability for rare 
realizations of the disorder which make H(uj) ~ and thus exp(— f3H) ~ 1 behaves like 
~ exp(— u 3 fi 2 /2g), for values of uj 3 » 2g/fi 2 (I have transformed their notation to ours). 
This could lead to the conclusion that the spatial probability distribution should also behave 
like a cubic power of to. This is true only if the probability distribution is completely 
dominated by the rare realizations of this kind, a fact that was never claimed by Villain 
et al.. It is certainly possible that other realizations which give lower values to exp(—/3H) 
but are nonetheless more abundant dominate the average value. Villain's argument is valid 
for any temperature. Our numerical results show that one certainly does not get a cubic 
behavior of the probability distribution at high temperature for values of uj » (g/jj 2 ) 1 ^ 3 . 
In that case the rare realizations of of the type considered by Villain still exist but they do 
not dominate the probability distribution. 

This is actually easy to understand. Let us define two different length scales in the 
problem. The first, introduced by Villain et al. which we call £1 is defined as 



with some undetermined constant A. The second length, which I introduced in eq.( |3.6| ). 




(5.2) 



is the length above which the probability for rare events of magnitude 1 goes like 



Aexp( 



(5.3) 




(5.4) 
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is the length for which the asymptotic behavior in the tail starts according to our investigation 
in Sec. 3. For to above this value we have found the behavior 

P(u) ~ exp(-/3/i^ 2 /2 + g p 2 \u\ + C{(3)) , (5.5) 

where C({3) is given in eq. (|3.8|). Since 

6 - 6 x ^ (5.6) 

We see that for T > 1.8, the cubic behavior is never realized since the contribution of the 
rare realizations of the Villain type to the averaged probability is smaller than the behavior 
given by eq.(5T5). This explain the perfect fit of the asymptotic behavior given in eq.( |4.3j ) 



and the numerical data, starting at a very low value of uj. 

On the other hand for T « 1.8 it is conceivable to have two tail regimes, the first with 
£i < uj < £2 in which the probability has the cubic behavior in uj because it is dominated 
by rare configurations of the Villain type, and a second tail regime for u > £ 2 for which the 
asymptotic behavior is changed to a quadratic behavior given by eq. (|5.5| ), if the results of the 
RS solution are valid (or possibly to a different asymptotic behavior which will emerge from 
a RSB solution). This is because one can easily check that at the border between these two 
regimes the two expressions given by eq.( |5.3| ) and eq.( [5.5| ) become comparable in magnitude, 
and the quadratic behavior wins for u > £2- As the (reduced) temperature is lowered from 
1 the first regime (cubic behavior) is expected to grow in size and include all of the tail at 
T = 0. Simulations of directed polymers at zero temperature, which can be mapped into 
corresponding results for the toy-model, |J showed the beginning of a change in the form of 
the spatial probability distribution from ~ exp(— ecu 2 ) to ~ exp(— cuj 2,2 ) at the onset of the 
tail, but did not go far enough into the tail region to confirm a cubic dependence. See also 



ref. [13]. We have previously observed cubic dependence in the tail of the directed polymers' 



spatial probability distribution at finite temperature. |15 
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There is some difficulty though, to explain the apparent cubic behavior in Fig. 3 purely 
in terms of Villain's configurations because we simulated "only" over 10 6 realizations and 
thus if the average is dominated by a single realization whose contribution is unity, the value 
of the average should be of the order of 10 -6 ~ exp(— 14) which far exceeds the value of 
the distribution in most of this region (except at the very beginning). One might argue 
that there are other realizations with somewhat smaller contributions to the average than 
Villain's which also gives rise to a similar cubic behavior but this needs to be verified. 

We hope that this work will stimulate further investigation of the behavior of the tail 
of the probability distribution at low temperatures. Three important questions that need a 
definite answer are 

1. Is there a RSB solution to our new variational equations? 

2. Is a cubic dependence of the log of the probability indeed realized over a large region 
when the temperature is very low? 



3. Is the asymptotic behavior given by eq. (|3.8| ) which works so well for high temperatures 
when uj > also valid asymptotically at low temperatures? 

We also hope that it will be possible to extend the variational method developed in this 
paper for the zero-dimensional manifold directly to higher dimensional manifolds in random 
media. 
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APPENDIX A: 



In this Appendix we derive eq. (|2.24 ) starting from eq. (|2.23|) . Let us shift the variables 
u a in eq.( |2.23| ) 

u a ->u a + \ a (Al) 
such that the linear term in u\ in the exponential is eliminated. This is achieved by choosing 

A a = -{i/(3)G la k (A2) 

We then find: 

Tr[(H - h)e~ (3h ] = J ^e ik " e"^ Gllfc2 jdujf- du n 

\ ab ab / 

6XP \ 2 l ^ abUJaUJ b] ■ ( A3 ) 

\ ab / 

We now expand the function f in a power series about 0, (this is used only as a tool derive 
our result. It may not be necessary for f to be analytic about 0) 

oo 

f(x) = J2fi xl (A4) 

1=0 

We use the following formula to integrate over u\ ■ ■ -uj n : 



x 



X 



J dux--- du n (co a - uj b ) 2s exp ( ~ J2(G ^ab^a^b J 

\ ab / 



(27r//3)«/ 2 (Det G) 1/2 ^!^f(Gaa + G bb - 2G ab y , (A5) 
provided s is a non-negative integer. For s being half-integer the integral is 0. We obtain 

f E / - G lb ) 2 k 2 /{3\ (G aa + G bb - 2G ab )/(3 ) exp Y.^U^a) } 

ab \ cd / J 

-$ + ^? G - + ^r( G "-"? G, -)( 1 -S)- (A6) 



where 

1 oo I / \ ' — s 

f(x, y) = -7= £ /i £ 9(h *n* + 1/2) f 
v 11 1 =0 s =o \ z yj 

9(l,s) =EE <W)/V f f] 2 J - m . (A7) 



j=0 m=0 

Since 



j / \m 



f°° dp e~v 2 (p + = £ *(J, «)r(a + l/2)q^\ (A8) 
we see that f(x, y) can be expressed in the form 



l r 00 , „2 „ / / x 



= e"^^= /°° cip /(2yp) e-^^v 7 ^) p . (A9) 

^/7T J-oo 

Using this representation for f(x, y) the integral over k in eq . (|A6|) can be performed, and we 
obtain the desired result given in eq. (|2.24j) . 
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Figure Captions 

Figure 1: Plot of the log of spatial probability distribution vs. uo 2 for T = 5 (solid line). 
The dashed line is the result of the hamiltonian variational approximation and the diamonds 
the results of our new variational scheme. 

Figure 2: Plot of the log of the spatial probability distribution vs. uj 2 for T = 1. The 
wiggly solid curve represent the result of averaging over 500000 realizations. The light 
wiggly curve is for 50000 realizations. The dashed curve and the diamonds are explained in 
the caption of Fig. 1. The solid smooth curve result from the asymptotic formula, eq.( |3.S| ). 

Figure 3: Plot of the log of the spatial probability distribution vs. u 2 for T = 0.23. 
The wiggly solid curve represent the result of averaging over 10 6 realizations. The dashed 
and dashed-dotted lines represent the hamiltonian RS and RSB approximations respectively. 
The diamonds are the results of the new variational scheme. The open triangles represent 
the empirical cubic approximation. 
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